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I. INTRODUCTION AND OVERVIEW 


Radio telecommunications, and the satellite services in particular, 
are governed by international agreements. The technical analyses on 
which such agreements may be based are performed by the International 
Consultative Committee on Radio (CCIR), using as input the work of study 
groups established for that purpose. NASA is actively supporting the 
work of these study groups in regard to satellite communications 
systems. The purpose of this grant is to develop computational methods 
and to perform engineering calculations which will assist in this task. 
The grant was initiated in January 1981 and augmented in July, 1981; 
this is the first summary report. 

Two specific tasks were undertaken during the period. The first 
deals with frequency-sharing between the Inter-satellite Service and the 
Broadcasting-satellite Service in the band near 23 GHz. The 1979 World 
Administrative Radio Conference (WARC-79) promulgated an Inter-satellite 
Services (ISS) band from 22.55 to 23.55 GHz. This band is to be shared 
with other services, including the Broadcasting-satellite Service (BSS) 
in the 22.55 to 23.00 GHz band [1]. Because of potential interference, 
the detailed implementation was not spelled out by WARC-79. Under the 
Grant, this Laboratory has assisted NASA in evaluating the potential 
interference between the BSS and ISS users, and in defining usable orbit 
geometries. The results will be found in Section II, which is intended 
as a preliminary draft for a paper to be submitted to an appropriate 
technical journal. Additional information, including computer program 
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documentation is contained in a technical report prepared under this 
Grant [2]. In its initial stages, this work also made significant 
contributions to a CCIR Study Group Report [3]. 

Since July 1981, the problem of synthesizing optimal and near- 
optimal spectrum and orbit assignments for broadcast satellites has been 
a second topic of investigation. The goal of broadcasting-satellite 
service system synthesis is to specify for each service area under 
consideration a set of channel assignments, polarizations, and an 
orbital slot in a manner which minimizes the amount of bandwidth 
required. Implicit in this is the requirement that a specified number 
of channels be supplied to each service area and that these channels be 
protected from interference. Protection ratios of 35 dB single-entry 
cochannel, 19 dB single-entry adjacent channel, 30 dB aggregate are 
typical. Other constraints, such as eclipse protection and minimum 
elevation angle, limit the flexibility of assignments. 

Our aim is to devise computerized spectrum/orbit synthesis 
techniques which will be useful at two important international 
conferences. The first, to be held in 1983, deals with broadcast 
satellite services (BSS), the distribution of information (especially 
television) to many users simultaneously, either as individual users or 
on a community or regional basis. The second, to be held in 1985, deals 
with fixed satellite services (FSS), i.e., point-to-point communications 
by means of satellites between specific fixed stations on the Earth. 
Since in both the FSS and BSS cases up- and down-link calculations are 
involved, the techniques to be developed for the 1983 and 1985 
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conferences exhibit a certain amount of commonality. However, there are 
also significant differences such as the number of up- and down-links to 
be considered, the coverage area (which affects the antenna patterns 
strongly), and the modulation methods to be employed. 

The spectrum/orbit assignment problem is still far from a 
definitive solution, but some progress has been made and is reported in 
Section III. It is our aim to have significant improvements over 
present techniques available for the 1983 conference, and a quasi- 
optimal computer code for the 1985 conference. 

Conclusions and recommendations are summarized in Section IV. 

II. SINGLE-ENTRY INTERFERENCE BETWEEN BROADCAST SATELLITE AND 

INTERSATELLITE SYSTEMS SHARING FREQUENCIES NEAR 23 GHz 

A. INTRODUCTION 

The 1979 World Administrative Radio Conference (WARC-79) authorized 
the use of the 22.55-23.55 GHz band by the Intersatellite Service (ISS) 
and the use of the 22.55-23,0 GHz band by the Broadcasting-Satellite 
Service (BSS) in region II, the Western Hemisphere [1], Full use of 
these bands requires frequency-sharing between the two services in the 
22.55-23.0 GHz frequency range. The purpose of this paper is to define 
acceptable satellite orbital assignments for this purpose. Calculations 
for some specific systems have been made previously by CCIR study groups 
[3-5]. Since the systems are still being defined, the more general 
problem is addressed here. Only the co-polarized case is considered; 
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additional discrimination is, of course, possible by the use of 

orthogonal polarizations, and the extegsiop to this case is 

' ■ < 

straightforward. Also we do not consider isolation which might be 
obtained by special modulation schemes, generally at the sacrifice of 
efficient use of the spectrum by at least one of the services. 

Under these conditions, the only available means for preventing 
interference between the two services is by the discrimination available 
from the antenna patterns and station geometries. Figure 1 shows the 
geometric parameters of interest, with the Earth radius exaggerated for 
clarity; the actual ratio between the Earth radius and the radius of the 
geosynchronous orbit is 1:6.6. The receiving (RX) ISS satellite is 
located at A, the transmitting (TX) ISS satellite is shown at C, the BSS 
satellite transmitter is at B, and a BSS Earth station receiver is shown 
at E. It is assumed that the ISS(TX) antenna is pointed at the ISS(RX) 
antenna, and vice versa, as required by good system design. Similarly, 
the Earth station antenna is assumed pointed at the BSS transmittejr. We 
also assume for simplicity that the BSS transmitter happens to be 
pointed directly at the Earth station for which interference is, being 
computed; while this need not be precisely true, it can be shown to 
have no significant effect on the conclusions to be drawn. 

In Figure 1, <|/i and t4 denote the angles between the axes of the 
receiving antennas and their respective potential interference sources; 
^ and describe the angles between the transmitter antenna axes and 
the receivers with which they might potentially interfere. Protection 
against interference results from the fact that, in general, receiving 
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Figure 1. Simplified geometry for analysis of ISS/BSS systems. 
(Earth radius exaggerated for clarity). 


antennas are not pointed at the interfering transmitters and 
transmitting antennas are not pointed at receivers where they would 
cause interference, i.e., the angles to 1^4 in Figure 1 are not too 
small. The task addressed here is to define "not too small" 
quantitatively in terms of the systems parameters and orbital 
assignments. 

Before looking at quantitative interference calculations, it is 
useful to consider the approximate geometries which make the various >I)-j 
angles small. The more complete coordinate system of Figure 2 is useful 
for this purpose. The symbols 0i, 63, 63 denote central equatorial 
angles with respect to the Earth. The center of the Earth is indicated 
by 0 ; the North Pole by P. Clockwise central angles are taken as 
positive, counterclockwise as negative, while with the angles only 
the magnitude Is of interest. In the figure, the ej central angle 
happens to be negative, and all other central angles positive, but this 
is not necessarily true. All the 0 angles can be defined over any 360 ° 
interval. The location of the BSS satellite is used as reference for 
all the central angles. The central angle from the BSS satellite to the 
ISS(RX) satellite is denoted by 0^, that from the BSS satellite to the 
ISS(TX) satellite by 63, and the equatorial angle from the BSS satellite 
to the Earth station longitude by 03. The Earth station latitude is 
denoted by a. The letters u,v,w,x,y,z denote distances. 

The specific equations relating these variables which are used to 
calculate the are given in the Appendix. From these equations, or 
alternatively from consideration of Figure 2 , the following conditions 
are evident. 
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To make ilg small, we require 61-82 near ± 160 ° and (03.il) near 
[(01+02)/2, 0 °]. The geometry for the + sign is shown in Figure 3 . 

In practice, the Earth station coordinates cannot approach (± 90 °, 0 °) too 
closely because of the requirement of a minimum elevation angle for the 
BSS as seen from the Earth station receiver. 

To make i|^ small, we require 0i near ± 160 ° and (03, il) near 
(± 80 °, 0 °) with the signs coordinated to be both positive or both 
negative. The corresponding geometry, for negative signs, is shown in 
Figure 4 . 

To make and tl;4 small, the requirement is that 62 be very small 
in magnitude. It should be noted that ilq decreases more rapidly than ^4 
as the magnitude of 02 approaches zero, as shown in Figure 5 . 

A condition which would result in maximum interference to the ISS 
system exists when and \p2 are both small simultaneously. This 
requires 02 small in absolute value and 0i near ± 160 °, and consequently 
I 02-01 1 will also be near 160 °. The geometry is shown in Figure 6; it 
is seen that all four of the angles become small under this condition, 
which represents a "pathological" situation where all possible 
interference contributions become large. In order to avoid its 
occurrence, it may be wise not to allow very long frequency-shared 
inter-satellite paths, e.g., paths longer than about 120° might be 
restricted to the part of the spectrum not shared with the broadcast- 
satellite service. 
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Figure 5. Geometry for small 



8. TYPICAL SYSTEM PARAMETERS 


The aim of this paper ifs to characterize acceptable geometries, 
from the point of view of mutual interference, for a considerable range 
of systems parameters, such as transmitted power, antenna gains, 
receiver noise temperature, etc. However, if the parameters are changed 
too drastically, the very nature of the problem changes. The "typical" 
system values about which variations might be made are shown in Table 1, 
together with an explanation of the symbols which will be used 
throughout this paper. The values for the BSS system are taken from 
CCIR Report 215-4[6,7], those for the ISS system from a proposed design 
[3]. The antenna discrimination of the ISS system is adopted from CCIR 
Report 558-1 [8] and is shown in Figure 7; those of the BSS system are 
taken from CCIR Report 810 [9] for community reception and 
WARC-BS-77-Annex 8 [10] for transmission. They are given in Figures 8 
and 9, respectively. Pointing and station-keeping errors are ignored in 
this study. Halfpower beamwidths are calculated from gain by the 
relationship 

4^ = /27 , 000/G (degrees) . (1) 

It should be noted that over most of their angular range, at angles 
well removed from the main lobe, the discrimination patterns are 
constant. This property turns out very useful because the number of 
variables is reduced; this facilitates the construction of families of 
universal curves for displaying the results in a particularly useful 
form. 
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TABLE 1 


TYPICAL SYSTEMS PARAMETER VALUES 

f - frequency - 22.75 GHz 

pT - BSS transmitter power - 34 dBW 
*^BSS 

- BSS transmit antenna gain - 36 dB 

‘"BSS 

- BSS receive antenna gain - 43 dB 

‘"BSS 

^BSS “ receiver noise temperature - 1100 K. 

^BSS ■ bandwidth - 40 MHz. 

^ISS " transmitter power - 10 dBW. 

_ iss transmit antenna gain - 52 dB. 

‘"ISS 

- ISS receive antenna gain - 52 dB. 

‘"ISS 

^ISS " receiver noise temperature - 1000 K. 

®ISS - ISS bandwidth - 125 MHz 
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Relative antenna gain (dB) 
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Figure 7. ISS antenna discrimination patterns. From 
references [8,10]. 
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Co'pohr and cross-polar reference patterns for receiving antenna 


Curve A: Co>polar component for individual reception without side-Iobc suppression 


0 

for 0 < (p < 0.25 <p. 

/(p \* 
- 12(--- 
\<P*/ 

for 0.25 <(). < <p < 0.707 <p„ 

- [ 9 . 0 +20 logio Q] 

for 0.707 ip. < q) < 1 .26 <p. 

-[8.5 + 25 logio (;-)] 

for 1 .26 <p, < <p < 9.55 (p„ 

- 33 

for 9.55 <p« < q> 


Curve A': Co*polar component for community reception without side>)obe suppression 
0 for 0 < <p < 0.25 

— 12 ^ for 0.25 tpo < (p < 0.86 tp» 


- [.0.5 + 25 log., (J) ] 


for 0.86 <p^ < <p up to intersection with Curve C (then Curve C) 


Curve B: Cross«polar component for both types of reception 
— 25 for 0 < (p < 0.25 (p* 

■ — ^30 + 40 logio j ^ — I I ^ for 0.25 <p« < 9 ^ 0.44 


- 20 

(30 + 2 ! logio 




for 0.44 <po < <p < 1.4 <p* 
for I.4tp* < <p < 2<p« 


— 30 until intersection with co-polar component curve; then as for co-polar component 


Curve C: Minus the on-axis gain 

Figure 8. BSS earth station (receive) antenna discrimination 
pattern. From references [9,10]. Curve A' is used 
in the calculations of this report. 
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Relative antenna gain (dB) 



0,2 0,3 0,5 


Relative 


Curve A: Co-polar component 


- 12 P- 

Wo 


-17.5 + 25 logio 


for 0 < <p < 1.58 q>o 
for 1 .58 (po < tp ^ 3.16 tpo 


for 3.16 tpo < <p 


after intersection with curve C: as curve C 


Curve B: Cross-polar component 

- ^40 + 40 logio ^ - 1 ^ for 0 <p 0.33 tp„ 


for 0.33 tp„ < tp 1 .67 tp„ 


- 33 


- ^40 + 40 logio I ‘ I j 

after intersection with curve C; as curve C 



Curve C: Minus the on-axis gain. 

Figure 9. 8SS satellite (transmit) antenna discrimination 
pattern. From reference. [10]. 
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C. INTERFERENCE CRITERIA 


Two commonly used criteria for interference calculations are the 
i nterference-to-noi se ratio I/N and the carrier-to-interference ratio 
C/I. The interference-to-noise ratio is useful in predicting the 
degradation of system performance due to interference directly by the 
relation 


C 

N+I 


C 

IT 


(1 




-1 


( 2 ) 


where C denotes carrier power, I interference power, and N receiver 
noise power. 

The ratio I/N is also useful to describe system acceptability in terms 
of an allowed noise margin N^i by the inequality 

N„ > 1 + X • (3) 

m 

On the assumption that the number of interferers to any one receiver 
will not be large, systems might be designed for I/N ratios in the range 
-5 dB to -10 dB. 

The carrier-to-interference criterion implies the notion that 
excess interference can be overcome by boosting carrier power, although 
at the penalty of possibly increased interference to other users. 

System performance can be calculated from C/I by 


= (C)-i + 



(4) 
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In terms of noise margin, the acceptability relation can be shown to be 
N (N -1) > . ■' (5) 

A more complete discussion will be found in the thesis by Wang [11], 


0. INTERFERENCE TO THE ISS SYSTEM 


1. The General Case 

Equations for the single-entry interference-to-noise ratio and 
carrier-to-interference ratio can be obtained by application of the 
Friis transmission equation [12,13] as 


I = 
N 


BSS 


^BSS '■’bSS 




BSS’ 2^ ^ISS ^^ISS 


>-^ISS’ l^c 


f2 (4 tt)2 x 2 k Tiss ^SS 


T 

C = PiSS Giss . 

I P gT 
BSS BSS 



x2 

z2 


( 6 ) 

(7) 


Since the goal here is to present the results in a form which is 
helpful for system design and specifically, wherever possible, by the 
use of universal curves, it is useful to eliminate as many variables as 
possible by combining them into universal factors which are easily 
calculated from the given system parameters. We define two factors Rj 
and \<2 


I . f2 . Biss 7i$$ _ ^ ^ c2 

^ PbSS ^ (4tt) 2 k 


^RSS'^BSS 






R .R ^gR 


BSS’ 2^ ISS ISS^'ISS 


'I' 1 
’ 1^ 


( 8 ) 
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(9) 


r GrCC _ r> 1 y2 

^ ' ’ “^ss(«^ss-*2)“5ss(^5ss>»i) ' ^ 

where the dependence on central angles arises through <li{6i,02)» 

x(®l) and z(0i,02). In application, the left sides of these 
equations, which define the universal factors R^, R 2 , are evaluated by 
the user from the system parameters, while the right side has been 
evaluated by computer and is presented in the form of graphs. In order 
to minimize the hand-computation labor, the numerical constants have 
been included on the right (computer-generated) side; in the case of the 
I/N calculation by Equation (8) this makes the equation dependent on the 
units employed, and the units Hertz, Watts, Kelvins are implied. While 
there are still too many variables involved in these equations to permit 
plotting universal curves, it is useful at this time to look at the 
shape of the allowed regions and relate them to the antenna 
discrimination angle conditions discussed previously. A computer 
code implementing the right side of Equation (8) has been written [14] 
and a plot appears in Figure 10. The forbidden region of unacceptable 
interference appears near, and is symmetric with respect to, the 
diagonal line 02 = 0, which slants from the lower right to the upper 
left of the figure. This is the condition for ti = 0. The broadening 
of the unallowed region near the center of the figure occurs because x, 
the separation distance between the BSS and ISS(RX) satellites becomes 
small so that unacceptable interference can be received even via the ISS 
receiving antenna si delobes. The broadening of the region near the top 
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PATHOLOGICAL REGION 


( ^3, L) ={ 60°,0®) 
R. = 254dB I/N = 


155 (T) -ISS (R) SEP 

( THETA 2 -THETA I ) 

Figure 10. Sample solution of interference to an ISS system 
by a BSS system, using the unrestricted computer 
code IING [14], I/N is used as criterion. The 
narrow regions indicate unacceptable geometries. 
THETA 1, THETA 2 refers to Q\,02 Figure 2. 

See Table 1 for system parameter values. 




left ( 02 - 0 i=- 16 O°, 01 = 160 °) occurs due to also being small, so that 
the BSS transmitter points at the ISS receiver. The combination of the 
two effects, >i small, small represents the pathological condition 
alluded to previously, which can be avoided by restricting the maximum 
length of f requency-shared intersatellite links to central equatorial 
angles on the order of 120°. 

The computer code used for producing Figure 10 is useful not only 
to give an intuitive feel for system behavior, but also as a definitive 
tool for evaluating allowable system geometries for proposed specific 
systems. It deals with all possible geometries. A similar presentation 
for the C/I criterion, based on Equation (9), appears in Figure 11. 

This computer code is also documented in the technical report [15]. 

2. Universal Curves 

When the ISS receiver is not illuminated by the main or near side 
lobes of the BSS transmitter, e.g., for the geometry of Figure 12, 
Equations (8) and (9) can be simplified by setting 

Sss ‘ ‘^BSS ^^BSS*'*'2^®1’®3’^^^ = 1 . (10) 

The reduction in the number of variables then allows the plotting of 
families of universal curves with the universal factors Rj or as 
parameter. One such family is required for each value of ISS receiver 
antenna gain. Figure 13 shows an example of such a universal curve for 
the case of I/N used as the criterion; note that for the systems of 
Table 1 with I/N = -10 dB, R^ has the value 274.1 dB. Figure 14 shows 
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Figure 11. Sample solution of interference to an ISS system 
by a BSS system, using the unrestricted computer 
code ICIG [15], C/I is used as criterion. See 
Table 1 for system parameter values. 
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curves of Figures 13, 14, 15 are applicable 
• (4»2 large). 


ISS (R) -BSS (T] SEP ( THETA 1 ) 

60. 120. 180. 240. 300. 3B0. 



( THETA2 - THETA I ) 

Figure 13. Universal cu^ve for interference to an ISS system 

by a BSS system when the ISS receiver is illuminated 
by the back or far side lobes of thd BSS transmitter. 
I /N criterion. 
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ISS (RJ -BSS (T) SEP ( theta i ) 

20.0 30.0 
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an enlargement of the most significant region of this plot (with the 
abscissa reversed). The set of curves corresponding to Figure 12 but 
with C/I used as criterion is shown in Figure 15. For the system of 
Table 1 with C/I = 30 dB, the value of R2 is 38 dB. The region of 
applicability of the universal curves of Figures 13, 14, and 15, as 
defined by Equation (10), is shown in Figure 16. 

In principle, universal curves, can also be constructed for the case 
when the ISS(RX) antenna receives interference only via its far side 
lobes, as illustrated in Figure 17, i.e., when is large and ^2 is 
small. Then analogously to Equation (10) we have 

*"ISS '^ISS ^^ISS’'^1^®2^^ = 0.1 , (11) 

which may be used to simplify the right side of Equations (8) and (9) to 
involve only the central angles, Gj, 63, and £. A computer code for 
this geometry for the I/N criterion has been written, but it turns out 
the there are no forbidden regions unless the value of R^ is on the 
order of 230 dB, some 40 dB below the values associated with the typical 
system of Table 1 [16]. The physical reasons are a) the ISS(RX) antenna 
gain typically is much higher than that of the BSS(TX) antenna, which 
must illuminate a substantial Earth region for broadcast purposes, b) 
the far side-lobe region for the BSS(TX) antenna is taken as 0 dBi 
[6,7], while that for the ISS(RX) antenna is taken as -10 dBi [8], c) 
for the interference geometry of Figure 12 the distance between the 
interfering transmitter at B and the affected receiver at A can become 
small, but this is not possible for the geometry of Figure 17, and 
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1 5S (R] -BSS (T) SEP { THETA i ) 

•SO. 120. 180. 2140. 300. 



Figure 15. Universal curve for interference to an ISS system by 
a BSS system when the ISS receiver is illuminated 
by the back or far side lobes of the BSS 
transmitter. C/I criterion. 
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THETR 1 (MRXl 

140 . 160 . 180 . 20 



( 20 “, 0 ® ) 
( 15®, 30® ) 
'( 1 5 ®, 0 ® ) 



( 0,,L) = 


( 15®, 30®) 
( I 5®, 0® ) 
( 20 ®, 0 ® ) 


B5S(T) GRIN 

Figure 16. Regions of validity for the universal curves of 
Figures 13 to 15. The universal curves are not 
valid for points between the respective contours 
showing earth station locations as parameter; 
however, for the region between the dashed lines 
the BSS transmitter is screened from the ISS 
receiver by > the Earth,^ so that interference is 
impossible. 



ISS(RX) 



30 


Figure 17. Example of a second geometry for which universal 
curves may be constructed (ij)i large). 






d) the smallness of is unrestricted, while that of is restricted 
by the requirement of a minimum elevation of the BSS transmitter at R as 
seen from the Earth station at E. Thus the geometry of Figure 17 does 
not lead to interference for systems parameter values even moderately 
close to those currently considered, and the curves are not needed. The 
same is true for the C/I criterion. 

E. INTERFERENCE TO THE BSS SYSTEM 

In analogy with Equations (8) and (9), the following relations can 
be obtained for the interference criteria when the ISS transmission 
causes interference in the Earth-station BSS receiver 

I f2 Tg Be ■ p ^ISS ojss (4s$»^)Ge i^E(fiE»M) , . (12) 

N Piss ^ (4Tr)2|< w2 

T 

c Piss Giss = r = i ^ , (i3) 

^ ’’bSS ^BSS ^ISS ^^ISS’’^3^’^E^^E*’*'4^ 

where the central -angle dependence occurs through 'PgfSi, 02. 03. 

'^4(02. 93.^). v(03,Jl), and w(02,03,Jl). 

Computer codes have been written to implement these equations, and the 
results for an example are shown in Figures IB and 19. respectively 
[17]. In these graphs, the region between the vertical dashed lines is 
excluded because the Earth would block the ISS link; the region above 
the top dashed line and that below the bottom one are excluded because 
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• ISS (T) -BSS IT) SEP 1THETA2) ' 

120.- -60. 0. BO. . 120. , . 180. 



I 55 (T) - I SS (R) SEP ( THETA 2 -THETA I ) 


Figure 18. Sample solution of interference to a BSS earth 
station by an ISS system, using the unrestricted 
computer code EING [17]. I/N is used as 
criterion. The small "spike" defines the 
prohibited region. The narrow vertical region 
between dashed lines defines geometries which 
are unrealistic because the ISS system is 
obstructed by the Earth. For regions above the 
top and below the bottom horizontal dashed lines 
interference cannot occur because the BSS earth 
station cannot "see" the ISS transmitter. 
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ISS m -BSS (T) SEP ( THETA 2 ) 

-120. -60. 0. 60. 120. IBO. 


®ISS “52dB 
^BSS = 

( 03 ,L) = (6O®,O«) 


J 


k R^ = 50dB 

m 


raO. 60. 120. 180. 240. 300. 3S0. 

I SS ( T) - I SS (R) SEP (THETA2-THETA I ) 


Figure 19. Sample solution of interference to a BSS earth 

station by an ISS system, using the unrestricted 
computer code EC I G [17]. C/I is used as criterion. 
See Figure 18 for details of the presentation. 
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the ;BSS transmitter would be below the horizon as seen from the Earth 
station. 

If the systems parameters for the typical system of Table 1 had 
been used for these calculations, no unacceptable region would have been 
shown in the graphs. To obtain the graphs, R 3 was lowered by 23 dB and 
R 4 was raised by 28 dB from their Table 1 values; this is equivalent to 
requiring I/N of -33 dB and C/I of 58 dB instead of the conventional -10 
dB, +30 dB values, with other parameters unchanged. Even so, the region 
of unacceptable interference lies in the quasi-pathological region; if 
satellite paths over 120 " in central angle length were assigned to the 
unshared portion of the spectrum, the unacceptable regions in these 
graphs would be covered by these exclusions. 

It should be noted that when 1^3 and i |;4 are both small in Equations 
(12) and (13), ij/i and ^ will also be small (see Section II, above), 
therefore a pathological geometry from the point of view of interference 
to the BSS system is also pathological with respect to the ISS system. 
From these calculations it appears that the interference to the ISS 
system is the more restrictive constraint, so that for a broad range of 
system parameters interference to the BSS system will automatically be 
within allowed bounds if pathological geometries are excluded and 
interference to the ISS system is kept within allowed bounds. Defining 
quantitatively the range of system variables for which this statement is 
true remains as a task still to be completed. 

When either or 1^4 is sufficiently large so that the interference 
involves only the far side-lobe region of the Earth station or the 
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ISS(TX) antenna, the problem can again be reduced to one amenable to 
universal curve presentation. Computer codes for these cases have been 
written and universal curves plots constructed [18], but the R 3 and R 4 
values had to be degraded even more unrealistically from the typical 
values of Table 1 before any unacceptable region appeared. 

F. CONCLUSIONS 

The most serious interference between proposed ISS and BSS systems 
near 23 GHz arises when the interference is transmitted from the main 
lobe or near sidelobes of the interfering transmitter and is received 
through the main lobe or near side lobes of the receiver. In this 
report this situation is termed "pathological". It can occur for both 
interference to the BSS system and to the ISS system. A necessary 
condition for this situation is a long ISS path. It is therefore 
recommended that long ISS paths (e.g., those subtending an equatorial 
central angle greater than 120'’) be restricted to the 23.00-23.55 GHz 
range, which is not frequency-shared with the BSS. 

When the very long ISS paths are excluded from frequency-shari ng, 
simple sets of universal curves can be constructed to define acceptable 
satellite and Earth station locations, based on either I/N or C/I as the 
acceptability criterion. To use these curves, the designer performs a 
simple multiplication of various system parameters to find a universal 
factor; contours of the universal factor then allow the range of 
acceptable Tocations to be read directly from the charts. 
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In principle, one set of charts is required to define interference 
to the ISS system, and a second set is required for interference to the 
BSS system; the acceptable region is then the intersection of the 
individually acceptable regions. In practice, interference to the ISS 
system turns out to be the more restrictive condition for a wide range 
of reasonable system parameters. In this case only the ISS charts need 
be consulted. The precise definition of the parameter range for which 
this statement is true is yet to be determined, but it appears to 
include most, if not all, practical systems. 

Programs have been written which allow a proposed solution to be 
tested by calculating the resulting interference criterion (I/N or C/I) 
explicitly. A proposed design procedure would include use of the charts 
to select a suitable geometry and then to verify it by the direct 
calculation of the resulting interference. 

III. BROADCASTING-SATELLITE SERVICE SPECTRUM/ORBIT ASSIGNMENT 
SYNTHESIS 

A. INTRODUCTION 

The goal of broadcasting-satellite service assignment synthesis is 
to specify for each service area under consideration a set of channel , 
polarization, and orbital slot assignments in a manner which minimizes 
the amount of bandwidth required. Implicit in the problem is the 
requirement that a specified number of channels be available to each 
service area and that these channels be protected from interference. 
Protection ratios of 35 dB single-entry co-channel, 19 dR single-entry 
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adjacent channel, and 30 d8 aggregate are typical. Other factors, such 
as eclipse protection and minimutn elevation angle, limit the flexibility 
of solutions. 

The BSS assignment synthesis problem is a discrete-continuous 
nonlinear optimization problem. The discrete design variables are the 
choice of channel assignments and polarizations; the continuous design 
variables are the orbital slot assignments. The problem is nonlinear in 
several respects, e.g., antenna patterns, angular calculations, and pre- 

/ : i 

detection interference-to-signal ratios. Finally, it is an optimization 
problem in the sense that the objective is to minimize the required 
bandwidth, subject to design restrictions such as eclipse protection, 
etc. 

Mathematically, this problem is extremely difficult if the stated 
objective is to find the optimum solution and supply proof of its 
optimality. Optimization theory can at most give an indication that a 
solution is a local optimum with respect to the continuous design 
variables. Some form of enumeration would be requi red to provide 
similar statements including the discrete design variables. Simpler but 
related combinatorial optimization problems that result from very strong 
assumptions are (in the language of computational complexity theory) NP- 
complete. This means they are among the hardest problems known. 

However, this pessimistic assessment refers to the possibility of 
obtaining a proven optimum. This is, of course, not necessary from a 
practical point of view. Indeed, the objective of this research is to 
develop methods that will give good (or acceptable) system synthesis 
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with a high degree of reliability and a minimum of computational 
expense. In this regard, efforts have been concentrated in four areas 
that seem to hold the keys to a successful approach. They are the 
following: 

1. Incorporation of the SOUP program, 

2. Development of user insights, 

3. Exact algorithms, 

4. Heuristic procedures. 

An illustration of how these four elements might be arranged in a 
synthesis process is given in Figure 20. It should be noted that we are 
still in the early stages of research in this area, and the approach may 
need to be changed as we progress. 

B. INCORPORATION OF THE SOUP PROGRAM 

The initial effort in this research program was to assess the 
appropriate role for the SOUP program in a synthesis method. It was 
assumed from the outset that SOUP would be the final judge as to the 
acceptability of a synthesis plan. However, it was not clear that an 
iterative method could afford to make extensive use of the entire SOUP 
calculation process due to limits on computation time. Therefore, after 
installing SOUP on our VAX 11/780 system, timing tests were run with 90 
test points. Run times on the order of 10 seconds confirmed our 
suspicions that any practical synthesis method could make only limited 
use of SOUP. 
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Figure 20. Synthesis process flow chart, 
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Attention then turned to examining the logical structure of the 
SOUP calculations in order to determine what subset could be included in 
a "stripped down" version to be employed to provide surrogate measures 
of the quality of the synthesis plan. This surrogate measure would be 
used to direct the iterative search procedures and therefore must be 
implemented as efficiently as possible. The heart of the SOUP 
calculation was found to be the calculation of certain antenna pointing 
separation angles. It proved possible to streamline the calculation of 
these angles significantly compared to their calculation in SOUP. 

For example. Figure 21 on the following page illustrates the 
calculation of the angle a between the vector from the satellite to the 
main (i.e., intended) receiver (M) and the vector from the satellite to 
a receiver suffering interference (I). If we let (J> refer to the 
satellite longitude, (|)i and 9i the longitude and latitude of I, (|>m and 
0M the longitude and latitude of M, and the orbital radius R = 6.6134 
earth radii then the modified calculation of a is as follows: 


lien = [ 2 {l-cos((|>j^-(|)j)cos( 0 ^^)cos( 9 j) - sin( 9 j^)sin( 9 j)}J^ ( 14 ) 

7 1 

ilail = [R + 1 - 2R cos(9|^)cos(iij^ - (15) 

p 1 

iibll = [R + 1 - 2R cos(9j)cos((t)j - (l))]^ (16) 

s = (lien + nan + iibii)/2 (17) 

a = 2 arcos[{s(s-neii)/( nan iibn)}?] . ( 18 ) 
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Figure 21. Illustration of angle calculation. 


41 


This calculation requires 14 function evaluations and 27 arithmetic 
operations as compared with 18 and 39, respectively, in the SOUP 
equivalent. Furthermore, it is possible to reduce these calculations 
even more by storing a few intermediate results, such as iieii, which do 
not change with satellite position. It should also be mentioned that 
this calculation deals with distances in multiples of earth radii while 
SOUP uses kilometers, which introduces a noticeable numerical error. 

Another aspect of the implementation of the SOUP program at OSU 
that required attention was discrepancies in the outputs for runs made 
with identical data at NASA Lewis Research Center and OSU. The 
differences showed up in the margin calculations with differences up to 
0.5 dB. In an effort to pinpoint the source of the error, runs were 
made on the VAX in both single-precision and double-precision modes and 
on the Amdahl 470 of the Instruction and Research Computing Center of 
OSU. None of the four sets of runs (NASA, VAX-single precision, VAX- 
double precision, and Amdahl) were in complete agreement. The tentative 
conclusion was reached that the differences are the result of 
trigonometric calculations with small angles. 

C. DEVELOPMENT OF USER INSIGHTS 

The second area of effort in this project was the development of 
user insights, i.e., a base of knowledge developed through experience 
with potential system designs. The rationale for this effort was 
twofold. First, to be able to design an iterative procedure for this 
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problem, one must possess some feel for how the system reacts to design 
variables, with respect to both direction and rate of change. Second, 

J 

due to the complexity of the problem it may not be possible to develop a 
method that can proceed with the synthesis of an acceptable plan without 
iteration with a user possessing substantial insight into the problem. 
For example, it is hard to imagine how a classical optimization approach 
could uncover the potential benefit of selective use of crossed-path 
geometry without user intervention. 

One of the approaches we have taken to develop this basis of 
experience is to implement a computer-assisted version of a manual 
synthesis method which is currently in use at NASA Lewis Research 
Center. At its present level of implementation, our program attempts to 
specify an assignment of satellite orbit locations that satisfies the 
discrimination criteria for every pair of service areas. The user 
inputs satellite orbit positions and the program calculates the 
discrimination between all pairs and provides this information in the 
form of a discrimination matrix. The user can then make changes in 
orbit positions and improve the matrix. Pairwise comparison is also 
provided as an option. 

The program does not, and is not expected to, provide optimum orbit 
assignments. However, it is intended to provide good starting 
solutions. It has already demonstrated the capability to provide 
important insights, such as highly unacceptable assignments and cliques 
of areas With little or no interaction. It also has provided a basis 



for testing the sensitivity of approximations in the angle calculations, 
which consume most of the computation time. 

An anticipated extension to the program is to include frequency 
assignments and polarizations and to reduce the demands on user input. 

D. EXACT ALGORITHMS 

1. General Considerations 

Although it is clear that no standard optimization algorithm can be 
applied directly to this problem in its entirety, standard techniques 
may have a role to play when some of the design variables have been 
fixed at trial values and the others are to be optimized. For example, 
suppose that each service area has been assigned a set of channels and 
polarizations and that now it is desired to determine the corresponding 
orbital slots that minimize total system interference. 

This problem can be attacked by classical gradient-search methods 
for continuous variable problems. This is, of course, precisely the 
problem form encountered under block allocation schemes. 

There is, however, a potentially large computational burden 
associated with this approach. A gradient-search algorithm requires 
estimates of the partial derivations of the objective function with 
respect to each of the design variables, in this case the orbit 
locations for each service area. Obtaining these estimates will require 
a number of system interference calculations equal to the number of 
service areas. Furthermore, due to the strong interactions that may be 
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anticipated between the variables, modified gradient procedures may be 
required to give reasonable rates of convergence. These will require an 
estimate of the Hessian matrix, i.e., the matrix of second partial 
derivatives, and the computations associated with this grow as the 
square of the number of service areas. Clearly, the use of SOUP with up 
to 90 service areas is out of the question since these calculations must 
be repeated for every trial solution. 

A promising approach to. overcoming the computational cost of a 
full-blown gradient search procedure is based on two considerations. 
First, it is reasonable to expect that most of the system interference 
will be isolated in relatively few, perhaps 7 to 1?, of the service 
areas. Therefore, only the orbit positions of these services areas and 

■ f 

those of the service areas causing the interference need to be included 
in the optimization. Second, this is another example where a surrogate 
measure of interference, one that can be calculated much more quickly 
than an exact calculation in SOUP, can be useful. 

As another example of the potential role of exact algorithms there 
are aspects of the overall problem that lead to combinatorial 
optimization problems. While it has not been determined precisely how 
the solutions of these problems fit into the total synthesis of the 
satellite communication system, it is clear that they may play on 
important part. 

Consider the situation in which each pair of regions has been 
specified as either interfering or non-i nterferi ng with respect to co- 
channel assignment. A pair of service areas are co-channel interferers 


45 



if the elliptical beam assigned to one service area reaches the other 
without sufficient angular antenna discrimination. This may be true in 
both directions, but need not be bidirectional for the pai r to be 
co-channel interferers. Clearly, service areas sharing a common border 
will be interferers, but geographic contiguity is not a necessary 
condition. 

We limit our discussion initially to the case of assigning a single 
channel to each service area so as to avoid co-channel interference. 
However, this is not as limited a scenario as one might suppose at 
first. The case of assigning a fixed number of channels to each service 
area does not change the nature of the problem since this may be viewed 
as assigning groups consisting of this fixed number of adjacent channels 
so as to avoid co-channel interference between groups. The consider- 
ation of adjacent -channel interference may be handled by selecting the 
channels for each group so that adjacent channel interference does not 
occur within a group. For example, group 1 could consist of channels 
1,3, 5, 7 and group 2 of channels 9,11,13,15 in the case of four channels 
per service area. Then groups 1 and 2 could be assigned to a pair of 
interfering service areas without causing adjacent -channel 
interference. 

The problem under discussion may be formulated as a graph-coloring 
in which the nodes of the graph represent service areas and a pair of 
nodes are connected by an arc if the corresponding pair of service areas 
are interferers. We begin by describing a set -covering approach 
suggested by Cameron [19] and then suggest alternatives. 
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2. Cameron's Approach 


Cameron considered the problem of assigning a single channel to 
each service area so as to minimize the number of channels required 
[19]. The restriction on the assignment pattern was that no pair of 
interfering service areas could be assigned the same channel, in order 
to prevent co-channel interference. He suggested solving a sequence of 
minimum-cardinality set -covering problems to determine the minimum 
number of channels required and the corresponding channel assignment 
pattern. 

Cameron observed that the problem could be formulated as a graph- 
coloring problem. However, since service areas need not be 
geographically contiguous to interfere, the graph need not be planar and 
hence the famous four-color theorem does not apply. Graph coloring is a 
notoriously hard combinatorial problem (even for planar graphs). In 
fact, finding a coloring that requires no more than twice the minimum 
number of colors is among the hardest combinatorial optimization 
problems. 

Although the minimum-cardinality set -covering problem used in 
Cameron's approach is also a very difficult combinatorial optimization 
problem, practical experience with problems of this type has been 
relatively good and several computer codes exist with promise of solving 
problems of the size encountered in the current context. 

To describe Cameron's formulation (which was not expressed 
mathematically), suppose we have a set of n service areas. We propose 
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using p channels for assignment to the service areas. We wish to 
determine 1) is p a sufficient number of channels to provide a feasible 
assignment? and 2) if so, what is the assignment ^pattern? 

We introduce E which is defined as the set of all interfering 
pairs of service areas, i.e.. 


E = {(q»r)l service areas q and r interfere}. 

Now define the following two sets of decision variables 

’ 1 if channel i is not used for service area j 

^ij ' 

0 otherwise 

’ 1 if channel i is used for service area j 

yij ' 

0 otherwise 


(19) 


(20) 

(21) 


for i=l, ..., p; j=l, ..., n. The channel assignment set-covering 
problem for p channels is then denoted CASC(P): 


p n P n 

minimize I I xji + I I yij (22) 

i=l j=l i=l j=l 


subject 

to 

e 

1 yij >1, j = 1. ...» n , 
i=l 

(23) 


Xiq + Xir > 1, i = 1, ..., p; (q,r)eE , 

(24) 


^ij + Yij > 1. i = 1» •«•» P; j=l. .... n . 

(25) 


^ij» ^ij = 0 Of* 1. 

(26) 
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The first set of constraints insures that each service area j has 
at least one channel assigned to it as defined by the y variables; 
the objective function will insure that the optimum solution will have 
exactly one channel assignment. The second set of constraints insures 
that, for each channel, any potentially interfering pairs of service 
areas must not assign that frequency to at least one of the two regions 
as defined by the x variables. The final set of constraints requires 
that for each (i,j) either the corresponding y variable or the 
corresponding x variable or both equal 1. Of course, if both equal one, 
they contradict each other and this is the essence of the determination 
of the sufficiency of p. 

Clearly from the second set of constraints, the objective function 
is at least np. If a solution to CASC(P) can be found with an objective 
function value equal to np, this demonstrates that p channels are 
sufficient since for any (i,j) 

’^ij ^ij ~ 

Thus the x and y variables are then consistent and the y variables 
specify a feasible channel assignment. However, if the optimum value of 
the objective function exceeds np then all feasible solutions have 

xij = yij = 1 7 (28) 

>T’- 

for at least one (i,j) and hence all solutions are contradictory. In 
such a case, one concludes p channels are not sufficient. 
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To determine the minimum nu’mbeh of channels required amounts to 
determining the smallest value of p for which CASC(P) has an optimum 
objective function value of np. Although the relative computability of 
set -covering problems is an attractive feature, there are some drawbacks 
to this approach.. These are the following: 

1. Several problems (for different values of p) may have to be 
solved. 

2. The problem is somewhat larger in the number of variables and 
in the number of constraints than in other possible 
formulations. 

3. In practice, many codes stop short of finding an optimum 
solution but can deliver good approximately optimal solutions. 
However, they are of no use in this approach. 

3. 0-1 Programming Formulation 


Assume as before that there are n service areas. Let m be a number 
of channels known to be sufficient. Trivially m=n will work; however, 
smaller values may be obtained by heuristically generating a reasonable 
(although probably not optimal) channel assignment. For any region j 
define Cj as the set of service areas that interfere with service area j 
and decision variables 


xi 


1 if channel i is used 
0 otherwise 


(29) 
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1 if channel i is assigned to service area j 

(30) 

(31) 

subject m 

to I yij =1, j = 1, n , (32) 

i=l 

m 

I Yij - nxi < 0, i = 1, .... tn , , (33) 

I yik ■'•.ICjl Yij < |Cj|, i=l, .... m; j=l, n, 

KeC-i 

(34) 

where 

Xi, yij = 0 or 1 , (35) 

and jCjl represents the number of service areas in the set Cj. 

The first set of constraints insures that for each service area 
exactly one channel is assigned. The second set of constraints insures 
that for each channel, if that channel is assigned to one or more 
regions, this fact is reflected by the corresponding x variable being 
equal to one. Finally, the last set of constraints insures, for each v 
channel and service area, that assigning the channel to a given service 
area eliminates the possibility of assigning it also to any interfering 
service areas. The objective function minimizes the number of channels 
used. 


[ 0 otherwise ' . i 
The problem can now be stated as 
m 

minimize I x^ 
i=l 
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Note that any solution to this problem, whether optimal or not, 
yields a feasible assignment of channels (unlike the Cameron formulation 
where for some p even optimal solutions may not yield a feasible 
assignment). Hence, approximately optimal solutions may be found, and 
these can be quite useful. Also observe the following comparison of 
problem size given in Table 2, where |E| is the number of pairs of 
interfering regions. Since m and p are essentially equivalent and |E| 
is potentially very large, the 0-1 program is considerably smaller then 
Cameron's. 


number of variables 
number of constraints 


Table 2 

Problem Size Comparison 
Cameron 0-1 program 


2 np 

nm + m 

np + n +*'p|E| 

mn + m + n 


4. An Alternative Set -Covering Formulation 


Consider again a problem with n service areas. Generate a set of 

subsets Pi Pq of these areas such that within any given subset no 

pair of service areas interfere. Insure that every region is contained 
in one or more of the subsets. Then define the parameters 
1 if service area j is in the subset Pj 
0 otherwise 


3ij - 


(36) 
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and the decision variable 


xi 


1 if subset P-j i's included’ in the cover 
0 otherwise 


(37) 


Consider now a solution to the following set -covering problem: 


minimize i (38) 

i=l 

subject q 

to I aij Xi > 1, j = 1, .... n , ■ (39) 

i=l ‘ ‘ ‘ ■ 

where 

x-j = 0 or 1 . (40) 

The solution selects a minimum number of subsets such that every 
service area is contained in at least one subset. This provides a 
channel assignment pattern in the sense that every subset selected for 
the cover is assigned its own channel. Service areas belonging to more 
than one subsets selected in the cover may be assigned a channel chosen 
arbitrarily from the several subsets indicated. 

The quality of the channel assignment (i.e., the number of channels 
used) depends, of course, on the selection of the set of subsets. If, 
however, the subsets are chosen as maximal -cardinality non-interfering 
subsets (i.e., no pair of service areas in the subset interfere and no 
further service areas can be added without destroying this property) and 
all such subsets are selected, the channel assignment given by the 
solution to the set -covering problem is guaranteed to use the minimum 
number of channels. 
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Unfortunately, the problem of determining all maximal -cardi nal ity 
non-i nterferi ng, subsets (maximal cardinality independent sets is the 
graph theory equivalent) is itself a difficult combinational problem. 

In practice, one would use a heuristic to generate a reasonable set of 
subsets in the hope of finding good and possibly optimum channel 
assi gnments. 

5. Multiple Channel Adjacent Frequency Model 

Although we believe that the previously considered single-channel, 
(co-channel) interference model addresses the multiple-channel adjacent- 
channel interference cases indirectly, we consider now a 0-1 model for 
addressing these issues directly. To this end we define Cj - the set of 
service areas that are co-^channel interferers to j, and C'j - the set of 
service areas that are adjacent-channel interferers of j. 

For simplicty of notation, assume that the service areas have been 
ordered so that the first n'<n service areas possess adjacent channel 
interferers. Also let rj be the number of channels to be assigned to 
service area j. Using the same set of variables as in the previous 0-1 
programming formulation the model is 

m _ 

minimize Z^x-j (41) 

i=l 

subject m 

to I yij = rj, j = 1, ..., n , (42) 

i =1 
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m 


(43) 


I yij - nxi <0, i = 1, 
j=l 

I Yik t ICjl Yii < |Cj|, i=l, .... m; j=l, .... n, 
KeCi 

(44) 

I (yi-lj + Yi+lj) + 2(C-j|yij < 2lC-jl 
KeC'j 

i = 1, .... m; j = 1 n' ; (45) 

i 

where 

j Yi j “ 1 • ' (46) 

The differences between this 0-1 model and the previous model are a 
more complex objective function and an additional set of constraints. 

The change in the objective function is required to insure that 
bandwidth, and not simply the number of channels used, is minimized. 

This was not necessary in the absence of adjacent-channel restrictions 
since whichever channels were used could be reordered to eliminate 
unused channels from the required bandwidth. This, of course, is 
generally not possible with adjacent-channel restrictions . Consequently 
the objective must insure explicitly that bandwidth is minimized. The 
additional set of constraints is needed to insure adjacent-channel 
restrictions are observed. 

It must be noted that these models are likely to be difficult to 
solve for a proven optimum solution for a 90 service area and 40 channel 
problem. However, obtaining good solutions may be relatively easy. 


55 



especially when only a subset of the service areas are considered, e.g., 
only service areas with unresolved interference problems. Our own 
computational experience and that of Baybars [20] confirms this. 

E. HEURISTIC PROCEDURES 

The review of three Canadian papers was an important part of the 
process of attempting to conceptualize the features of heuristic 
approaches likely to provide a good starting solution. 

Chovinard and Vachon [21] present a method based on exhaustive 
(but implicit) enumeration of all possible channel assignments and 
polarizations to service areas, given a preassignment of orbital 
position to service areas. This brute-force approach seems unlikely to 
be of much help in realistically sized problems and does not offer any 
insights into the problem. Nedzela and Sidney [22] offer an approach 
that is a heuristic based on matrices indicating the freedom of choice 
for the service areas remaining to be considered later when assigning a 
channel, polarization, and orbital slot for a service area under 
consideration currently. By making the current assignment on the basis 
of maximizing the resulting freedom for later choices, a sequence of 
assignments is made that either results in a successful plan based on a 
prespecified number of channels or a procedure failure. At each step, 
the selection of service area to be considered next is made by choosing 
the one with minimum remaining freedom. Christensen's procedure [23] is 
an interactive system that offers the user a menu of routines, some of 
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which are automatic and some of which require manual input, to aid in 
synthesizing a plan. Both this paper and the Nedzela and Sidney paper 
offer some wel 1 -concei ved heuristic approaches which are part of our 
current thinking on the subject. 

The drawbacks of the Canadian procedures are twofold: 1) some of 
thex approximations that have been made may not be valid. For example, 
in the method of Nedzela and Sidney, assignments based on the minimum 
freedom matrix may not be feasible when tested by SOUP and, of course, 
there is no guarantee of optimality; and 2) a combination of all three 
methods was proposed with a supervisory iterative process. However, 
there seem to be certain difficulties encountered in trying to combine 
the methods and, so far, integration of the three procedures has not 
been achieved. Consequently, though these procedures are certainly 
useful, they fall short of solving the assignment synthesis problem. 

IV. CONCLUSIONS AND RECOMMENDATIONS 

The problem of allocating satellite orbital locations to the 
Broadcasting-Satellite Service and the Inter-satellite Service near 23 
GHz has been solved with respect to single-entry interference, subject 
to the constraint that very long inter-satellite links (e.g., those 
separated by more than 120° of equatorial arc) be allocated frequencies 
in the unshared portion of the band. The allocation procedure involves 
the use of universal charts. To use these charts, the designer 
calculates a universal factor by simple multiplication of certain system 
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parameters; acceptable geometries are then found from contours of that 
parameter value on charts whose coordinates are satellite separation 
angles. Programs for checking the validity of the resulting assignments 
have also been prepared. 

It is recommended that an atlas of such charts be prepared and 
published, together with instructions for its use. Certain minor points 
remain still to be resolved, e.g., whether the assignment on the 
basis of interference to the ISS always assures protection also to the 
BSS for practical system parameter values, or whether two sets of charts 
may be required in certain cases. 

Considerable progress has been made in formalizing the concepts of 
broadcasting-satellite service assignment synthesis. An important 
feature of our suggested approach is the combination of exact 
algorithms, heuristic procedures and user insights. We are confident 
that the formulations and solution procedures described in this report 
are effective methods for the corresponding subproblems of the overall 
system synthesis problem. However, the most challenging part of the 
problem remains to be formally addressed. This is the design and 
testing of a software package that combines and effectively interfaces 
the several subprocedures already developed. This effort is ready to be 
initiated. In the process, additional procedures and approaches may 
emerge. 


58 




REFERENCES 


[1] Final Acts of the World Administrative Radio Conference (1R79) ^ 
Radio Regulations, Part A» Chapter III, Article ft. Section IV, 
p. 140 (RR8-156), International Telecommunication Union, Geneva, 

1979. 

[2] Wang, C.W. and Levis, C.A., Allowed Orbit Geometries for Frequency 
Sharing by Broadcasting Satellite and Inter-satellite Systems near 
23 GHz. Report 713533-1, The Ohio State University ElectroScience 
Laboratory, Department of Electrical Engineering; prepared under 
Grant No. NAG 3-159 for NASA-Lewis Research Center, 21000 
Brookpark Road, Cleveland, Ohio, February 1983. Also available, 
with some changes, as M.S. Thesis, The Ohio State University. 

[3] International Radio Consultative Committee, CCIR Study Groups, 
Period 1978-82, "Sharing Between the Inter-satellite Service and 
the Broadcasting-satellite Service in the Vicinity of 23 GHz". 
Proposed Modification to Draft Report AC/10-11, Study Programme 
2P/4. Document lO/llS/136 E (USA), dune 9, 1981. 

[4] International Radio Consultative Committee, CCIR Study Groups 
Period, 1978-1982, "Sharing Between the Inter-satellite Service 
and the Broadcasting-satellite Service in the Vicinity of 23 GHz". 
Draft New Report lO-llS-2, Document 10-11S/96-E, October IB, 

1980. 

[5] International Radio Consultative Committee, CCIR Study Groups, 
Period 1978-82, "Protection of the Broadcasting-satellite Service 
from Interference by the Intersatellite Service in the Vicinity of 
23 GHz". Document 10-11S/115-E (Japan). 

[6] International Radio Consultative Committee (CCIR) Report 215-4, 
"Systems for the Broadcasting Satellite Service (Sound and 
Television)". Recommendations and Reports of the CCIR, 1978, 
XlV-th Plenary Assembly, Kyoto, 1978, Vol . XI, pp. 163-185. 
Published by International Telecommunication Union, Geneva, 1978. 
ISBN 92-61-00761-9. 

[7] International Radio Consultative Committee, CCIR Study Groups, 
Period 1978-82, "Proposed Modifications to Report 215-4" 

(lO-llS-3) Document 10-11S/104-E, p. 13, column 10, October 14, 
1980. 


59 



[8] ; International Radio Consultative Committee (CCIR) Report 558-1, 

• "Satellite Antenna Patterns in the Fixed Satellite Service". 

Recommendations and Reports of the CCIR, 1978. XlV-th Plenary 
Assembly, Kyoto, 1978, Vol . IX, pp. 243-262. Published by 
International Telecommunication Union, Geneva, 1978. ISBN 92-61- 
00741-4. 

[9] International Radio Consultative Committee (CCIR) Report 810, 

"Broadcasting-satellite Service (Sound and Television): Reference 

Patterns and Technology for Transmitting and Receiving Antennae". 
Recommendations and Reports of the CCIR, 1978, XlV-th Plenary 
Assembly, Kyoto, 1978, Vol. XI, pp. 271-285. Published by 
International Telecommunication Union, Geneva, 1978. ISBN 
92-61-00761-9. 

[10] Final Acts of the World Administrative Radio Conference for the 
Planning of the Broadcasting-satellite Service in Frequency Bands 
11.7-12.2 GHz (in Regions 2 and 3) and 11.7-12.5 GHz (in Region 
1), Geneva, 1977. International Telecommunication Union, Geneva. 
ISBN 92-61-00491-1. 

[11] Wang, C.W. and Levis, C.A., op. cit. (Chapter II). 

[12] Wang, C.W. and Levis, C.A., op. cit. (Chapter IV). 

[13] Kraus, J.D. and Carver, K.R., Electromagnetics, 2nd ed., McGraw- 
HilU 1973, pp. 696-699. 

[14] Wang, C.W. and Levis, C.A., op. cit. (Appendix D). 

[15] Wang, C.W. and Levis, C.A., op. cit. (Appendix H). 

[16] Wang, C.W. and Levis, C.A., op. cit. (pp. 31-37 and Appendix G). 

[17] Wang, C.W. and Levis, C.A., op. cit. (Chapter IV, Section B and 
Appendices J and M). 

[18] Wang, C.W. and Levis, C.A., op. cit. (Figures 20, 21, 23 and 
Appendices K,L,N). 

[19] Cameron, S.H., "The Solution of the Graph-Coloring Problem as a 
Set -Covering Problem", IEEE Trans. Electromagnetic Compatibility, 
-Vol. EMC-19, No. 3, August 1977. 

[20] Baybars, I., "Optimal Assignment of Broadcasting Frequencies", 
European Journal of Operational Research, Vol. 9, No. 3, March 
1982. 


60 



i 

[21] Chouinard, G. and Vachon, M., "A Synthesis Procedure of a Plan by 
Computer Using an Assignment Method and Inductive Selection", ITU/ 
CITEL/DOC-Canada Preparatory Seminar for 1983 RARC, Ottawa, 

Canada, May 1981. 

[22] Nedzela, M. and Sidney, J.R., "A Synthesis Algorithm for 
Broadcasting Satellite Service", ITU/CITEL/DOC-Canada Preparatory 
Seminar for 1983 RARC, Ottawa, Canada, May 1981. 

[23] Christensen, vl., "BSS CAPS, A System Description", ITU/CITEL/DOC- 
Canada Preparatory Seminar for 1983 RARC, Ottawa, Canada, May 
1981. 


61 



APPENDIX 


EQUATIONS FOR CALCULATING DISTANCES ? 

AND ANTENNA ANGLES FROM CENTRAL ANGLES 

The equations below may be used to calculate the distances and 
angles shown in Figure 2 . The symbols $ and r denote the radius of the 
geostationary orbit and of the earth, respectively. 


u^ = 0^ + r^ - 2r0 cosjl cos(02-Q^) (A*l) 

= 0^ + r^ - 2r0 cos£ cosG^ (A*2) 

w^ = 0^ + r^ - 2r0 cos«- cos (02-62) (A. 3 ) 

X = 20 sin 1^1 (A. 4 ) 

y = 20 sin \^\ (A« 5 ) 

z = 20 sin (A. 6 ) 

2 

1^2 = cos“^[(x^+z^-y^)/( 2 xz)] (A« 7 ) 

1^2 = cos~^[(x^+v^-u^)/( 2 xv)] (A* 8 ) 

i|^2 = cos“^[(w^+z^-u^)/( 2 wz)] (A. 9 ) 
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(A. 10) 


= cos“H(w^+v^-y^)/(2wv)] 

, 0 > 0 (A. 11) 

1 2 12 

= 180° - „ 0^ < 0^ (A. 12) 

Note that in Equations (A. 11) and (A. 12) the ranges of 0i, 02 are 
from 0° to 360°. 
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